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Abstract: 

We compute the leading logarithmic behaviour of the cross-section for the production of a pseudoscalar 
Higgs boson in gluon-gluon fusion to all-orders in perturbation theory, in the limit of large partonic 
centre of mass energy. We also calculate the Higgs rapidity distribution to the same accuracy. We 
include the contributions of top and bottom quarks, together with their interference. Our results are 
given in terms of single and double integrals, evaluated explicitly up to next-to next-to leading order 
(NNLO). We use our results to improve the known NNLO inclusive cross-section computed in the 
effective theory where the fermions in the loop are integrated out. The size of finite fermion mass 
effects on the inclusive cross-section is found to be small, reaching a few percent only for large values 
of the pseudoscalar mass. 

1 Introduction 



The first collisions at the Large Hadron Collider (LHC) at CERN marked the beginning of a new era in 
particle physics. The four primary experiments at the LHC have started taking data at energies never 
explored before. Their main target is to understand electro-weak symmetry breaking and the origin 
of mass in the universe. Over forty years after its proposal, the Higgs mechanism for giving mass to 
fundamental particles will finally be tested. Supersymmetric theories generally predict a richer Higgs 
sector than the Standard Model (SM). Li the Minimal Supersymmetric Standard Model (MSSM) for 
instance one introduces two complex Higgs doublets, which originate five physical Higgs bosons: two 
neutral scalars (/i, H), two charged scalars {H^) and, finally, a neutral CP-odd state {A). 

The most important production channel for the pseudoscalar Higgs boson at hadron colliders is 
gluon-gluon fusion, via a fermion loop, and associated production with jets. However, the phenomenol- 
ogy is much richer than the case of the SM Higgs, because the coupling of the pseudoscalar to the 
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up-type quarks decreases with tan^ = vi/v2 (the ratio of the two vacuum expectation values), while 
the one to down-type quarks increases. As a consequence the coupling to the 6-quark becomes impor- 
tant and eventually dominates the top contribution for tan /? > 10 (for a recent review see for instance 
Ref. [1]). 

In collinear factorization, the production cross-section for gg ^ A can be written as 
a{as;Xh,mA;mt,mb) = /" dxi f f \ ^^^^^ ^-^^ 

. J Xi J X2 \ X1X2 J 

where x^ = m\/ s and both renormalization and factorization scales have been set equal to the pseu- 
doscalar mass m^. The partonic cross-section aij has been computed in QCD at next-to-leading order 
(NLO) [2]: 



fj'=t,b 



S,^Sil-x) + ^cff (x;^,^)+Oia^^) 
TT \ ruA ruA J 



(2) 



The NLO corrections turned out to be large (0(70%)) and motivated a NNLO calculation. However, 
a full treatment would require the evaluation of 2 ^ 1 diagrams at three loops with massive lines. For 
this reason the NNLO corrections have been computed in an effective theory [3] in which one integrates 
out the contribution of the heavy fermion running in the loop [4, 5, 6]. This approximation at NLO 
works very well, far beyond its naive range of validity, essentially because the cross-section is dominated 
by the emission of relatively soft gluons which cannot resolve the fermion loop. However, in the MSSM 
this approximation is more delicate because the mass of the pseudoscalar can be much larger than the 
top mass. Furthermore, at large tan /3, the coupling to the 6-quark becomes increasingly important and 
the use of the effective theory less justified. 

It is therefore important to quantitatively asses the role of finite fermion mass corrections beyond 
NLO. In the case of SM model Higgs production, recent analyses [7, 8, 9] confirm that also at NNLO 
finite fermion mass effects are negligible. Those analyses are based on an asymptotic expansion of the 
Feynman diagrams in inverse powers of the heavy fermion mass nif. However, it was soon realised [10] 
that this expansion fails to converge in the limit of large partonic centre of mass energy s (i.e. small 
x = JT^Higgg/s), and it gives rise to a spurious power growth. Fortunately the small-x bevahiour of the 
partonic coefficient functions can be computed using a technique called A;T-factorization [11, 12]. This 
has been done in the case of the SM Higgs in [13, 14] and combined than with the asymptotic expansion 
in Refs. [7, 8, 9]. 

In this letter we explicitly compute the small-x behaviour of the partonic cross-section for gg ^ A 
at NNLO QCD and we match it to the effective theory result, which can be viewed as the zeroth order 
in the asymptotic expansion. This enables us to give a first estimate of finite fermion mass effects in 
NNLO pseudoscalar production. 



2 The off-shell cross-section 

We use fcr-factorization [11, 12] to determine the leading logarithmic behaviour at small-a; (LLx) of 
the partonic coefficient functions to any desired order in the strong coupling. We start by considering 
gluon-initiated subprocesses only, as the other partonic channels can be derived using colour-charge 
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relations. Thus, we compute the leading order cross-section for the production of a pseudoscalar Higgs, 
keeping the incoming gluons off their mass-shell 



g*{k-i,n,A,e-i) + g*{k2,i',B,e2) — > A{q), (3) 

where A, B are colour indices and are the polarization vectors, which will be specified later on. The 
leading order amplitude is obtained by computing the one-loop massive three-point function 

M = Y,[MkiM) + Af{k2M)], (4) 
/ 

where the index / runs over the fermionic flavours in the loop. In our study we consider f = t,b. 
The coupling of the pseudoscalar Higgs is proportional to the mass of the fermion running in the loop: 
—i—{^75, where vq = \/vi +V2 and gAf is a dimensionless quantity, which in the MSSM, for instance, 
is gAt = l/tan/3 for up-type quarks and gAb = tan/3 for down-type ones. The calculation of the 
amplitude is then straightforward; in = 4 — 2e we have 



Vo 



d^l T r[(/ + m/) M/- M + mf)j5{/+ h + ruf) M] 
(l^ - m2.) (^(l - kiY - m^) {{I + k2Y - mfj 



The treatment of 75 and of the Levi-Civita tensor away from four dimensions is in principle a delicate 
problem [15, 16, 17]. This computation is analogous to calculation of the axial anomaly. It is well known 
that, although the result of the loop integral is finite, there are some issues regarding the regularization 
procedure one must take care of. With this in mind, using standard manipulation of Dirac matrices 
and Feynman parametrization of the loop integrand, we obtain: 

Af = -4^mfgAfe{kl,k2,el,e2){gs^i^^f—{TRS^^)x 

Vo 

^ ^dadf3 ' ^^^') _ 



^^^^ P - (mj - a(l - a)kl - /3(1 - /3)A;2 - 2a^ki ■ k-^ 



where 



e{ki, k2,ei, 62) = e^t^pak'^k^ele^ and rr = ^\^(/1^2g)"^ • C'') 

We notice that this result has a simpler form than the one for the production of a scalar Higgs (see for 
instance Ref.[13]) and it is basically given by the one-loop three-point function [18]. Since everything 
is regular, we can safely go back to D = 4 dimension and we obtain 

M=J2^^f=Il S^''^^<k^,k2,e^,e2)hiklklq^m},m},mj). (8) 
f=t,b f=t,b ^ '"^ 
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In order to determine the high-energy behaviour of the cross-section it is convenient to write the 
kinematics in terms of: 

ki = ZiPi + kj, = 0) ~ ^ = 1) 2, 

ml = q^ = {ki + k2f, s = 2pi-p2; (9) 

we also find useful to introduce the following dimensionless ratios: 

x = ^, rf = ^, ^i=^^,i = l,2. (10) 
s rrij^ 

In this notation the one-loop amplitude looks like 

M=y^ (5^^^e(A:i,fc2,ei,e2)47rr/Co(ei,6;T/), (11) 

where to make contact with the case of a scalar Higgs boson of Ref.[13], we have introduced the 
dimensionless form factor 



2 

m 



(12) 



In order to obtain the off-shell cross section in ^T-factorization, we specify the polarization vectors 
Cj to be 

we then square the matrix element Eq.(ll), sum (average) over final (initial) state colour and spin and 
multiply by the phase space factor dV, obtaining 



r(x,Ci, 6; rt,r6) = / dV\M? 

= 27r3a2GFV2^'"^{5AtV2|Co(6,e2;rt)|' + 5V2|Co(ei,6;T6) 

}j 
■ 

X5 - 1 - ^1 - ^2 - 2^/66 COS (^^ , 



sm^ 
Ax 



(14) 



where we have introduced the Fermi constant Gp = the azimuthal angle (p between the 

transverse components of the incoming momenta. The small-x leading logarithms of the collinearly 
f actor ized coefficient function can be computed by considering an impact factor, which is defined as the 
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triple Mellin transform of the off-shell cross-section Eq. (14): 

pl poo poo 

Jo Jo Jo 

+25At5A6rtr;.ReCo*(ei,6;rt)Co(6,6;r6)| T"^- ''''' Hv' (15) 

where we have performed the AT MeUin using the delta function. We remind the reader that the above 
impact factor leads to a result for the coUinear factorized coefficient functions which is not in MS. 

The mismatching between the two factorization schemes is known and described by an overall factor 
R [12]. However, we are interested in the expansion of the resummed result at NNLO and the scheme 
dependence starts one order higher: 

R=l + 0{al). (16) 

Therefore, at the accuracy we are working at, we can safely ignore all issues concerning the scheme 
dependence of our results. 



3 Determination of the high-energy behaviour 

The leading small-x singularities of the coefficient function can be now determined from the impact 
factor h{N,Mi,M2;Tt,Tj,), identifying Mi = M2 = 'ys{as/N), the anomalous dimension which is dual 
to the leading order BFKL kernel [19]: 

and then expanding in powers of ct^ (i.e. in powers of Mi, M2), taking the — t- limit. Before 
computing the M Mellin transforms in the full theory, it is instructive to look at the result in the heavy 
fermion approximation. In this effective theory we expect a spurious double logarithmic behaviour at 
small-x [11, 20]. We investigate this case in the following section. 

3.1 Heavy fermion approximation 

We consider the case where there is only one fermionic fiavour running in the loop, which is infinitely 
massive. We have 

^i^^CUi„h;r,) = -^ + o(±y (18) 
so the impact factor becomes 



d^2e2 



f^'^ dip sin^ (p 

Jo 



27r / 9 nnr- \ ^ ' ^^^"^ 
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Prom the expression above it is clear that one cannot set N = 0, because the integrals would then 
diverge for every value of Mi, M2. The situation is totally analogous to scalar Higgs production and 
double high-energy logarithms appear to all orders in perturbation theory [20]. The angular integration 
can be performed in terms of hyper geometric functions, but it is easy to see that it is actually safe to 
set iV = inside the angular integral. The azimuthal integration then trivially gives: 



poo f'OO 

h{N,Mi,M2;oo) = a^ d^i^^'-' d^, 

Jo Jo 



M2-1 



[l + OiN]], 



(1 + ^1 + 6)^ 

where we have introduced the leading order cross-section in the heavy fermion limit 



(20) 



-0^ = 



2567r 



(21) 



We note that the pseudoscalar cross-section differs from the scalar one by a numerical factor 9/4. This 
is expected as the two effective Lagrangians are 



eflf 



Stt 



(22) 



In the heavy fermion approximation, once we have factored out the leading order result, the high-energy 
leading (double) logarithms are the same for scalar and pseudoscalar Higgs: 



h{N, Ml, M2; 00) 



-0^ 



-0^ 



1 + 



M1 + M2 (Mi+M2\ 



N 



+ 



— ; 



+ 



a,2CU (2CaV 



The LLx result in x-space is easily obtained by inverting the Mellin transform: 
agg{x; 00) = [(5(1 - x) + ^2Ca In - + (^) ' ^ In^ - + 0(«f ) . 



(23) 



(24) 



So far we have been considering the gluon-gluon channel only. The generalization to all the other 
partonic subprocesses is straightforward in the high-energy limit [11, 12]. We obtain 



aqg{x;oo) = 

aqq{x;00) = Uq 



TT X V TT / 

\TrJ 6! 



6! 



'2Ci-ln3i+0(a3) 



(25) 



The above results Eq. (24), (25) agree with the small-x limit of the known NNLO coefficient functions [5]. 
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3.2 Finite fermion masses 

We now go back to Eq. (15) in order to compute the LLx behaviour in the full theory, with top and 
bottom quarks running in the loop. We first notice that the presence of the form factor |Cop ensures 
that the M Mellin transforms have finite radius of convergence when = 0. We write the impact 
factor as the sum of different contributions: the ones coming from the square of each diagram and then 
interference ones: 



h{0,Mi,M2;Tt,n)= H0,Mi,M2;Tf,Tf). 
f,f'=t,b 

where we have introduced: 

h{0,Mi,M2;Tf,Tf,) = 

77— . — ^ — M1M2 j dCiCi / "^2^2 <^o(a>?2; -7)00(^1, C2;'r/') 

\l + Offi] Jo Jo 



(26) 



2\/2 (1 



i + s 



ff 



POO J r^*i 

Jo ^0 



d Co*( 6,0;T/)Co(gi,0;T/. 
(O,O;r/)Co(O,O;r/0 



--{Ml + Ml) I d^ln^a 



2 . Co*(6,0;r/)Co(6,0;r/. 







Co*(O,O;r/)Co(O,O;T/0 

5' Co*(a,6;r/)Co(a,6;r/0 



f-OO f'OQ 

+M1M2 d^i d^2ln^i\n^2^. „ ^nm n \ 

Jo Jo d^id^2 Q(0,0;r/)Co(0,0;r//) 



+ ... 



with 



^0 iV^V) = ^ ^Co(0,0;r/)Co(0,0;T/0- 

In particular when r/ = Tf , the LO contribution to the cross-section is 

<^oirf,rf) = ^ ^ |go(0,0;r/)P = -^£-\rf f {rf)\' , 



2v/2 



(27) 



(28) 



(29) 



firf) 



In 



l-Tf 



ITT 



arcsm 



if Tf<l 

if Tf > 1. 



(30) 



The high-energy behaviour is then determined by the identification Mi = M2 = 7s {ag/N). The 
result in momentum space is then easily recovered by inversion of the AT-Mellin transform. We obtain 



Crgg{x;Tt,n) 



a 



S(l-x) + ^2CACliTf,Tf 



TT 



+ {^ycl (2c2(r/,r/0 + cn(r/,r/0) In^ + 0(«^) 



(31) 
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where the coefficients are expressed in terms of the following integrals 

d C*(6,0;r/)Co(a,0;r^O 



ci(r/,r//) 
C2(t/,t// 

Cu{Tf,Tf 



"^^'^""^'d^ Co*(0,0;T/)Co(0,0;r/0 ' 



1 TjA d C*{^uO;rf)Co{^i,0;rr) 

2 Jo ^'d6 Co*(0,0;r/)Co(0,0;r;0 ' 

'^Vo Co*(0,0;r,n(0,0;r,,) ' 



(32) 



(33) 



We evaluate numerically the above integrals; the results are tabulated in Table 1 and 2. We use 
771-5 = 4.19 GeV and mt = 172.0 GeV. The LLx behaviour of the other subprocesses is easily obtained 
with colour-charge relations: 



Crqg{x,Tt,Tb) = 



f,f'=t,b 



1 + Sff 



--^CFCi{Tf,Tfi) + i^-^j CaCf (c2(r/,T//) + cii(r/,r/')) In - 



= E ■ 

f,f'=t,b 



1 + Sff 



+ 



(34) 



We note that when the fermion mass is kept finite, the cross-section for pseudoscalar Higgs produc- 
tion has the expected single-logarithmic behaviour. The coefficients of the LLx singularities differ from 
the scalar Higgs case [9]. As a check of our procedure we have compared our asymptotic behaviour to 
the small-x limit of the full NLO calculation of Ref. [2] and we have found perfect agreement. As an 
example in Fig. 1 we plot the 0{as) coefficient function in the quark- gluon channel Cql^ in the case 
of just a top quark running in the loop, as a function of x. The green dashed curve is the effective 
theory result, which exhibits the spurious logarithmic growth at small x, while the solid curves are for 
mt = 172.0 GeV and niA = 100, 200, 300, 500 and 1000 GeV respectively. When the fermion mass in 
the loop is kept finite the coefficient function has the expected constant behaviour at small x, which 
is in numerical agreement with the coefficients in Table 1, once multiplied by the appropriate colour 
factor Cf- 



3.3 Rapidity distribution 

The formalism which enables one to resum high-energy enhanced contributions has been recently ex- 
tended to the case of rapidity distributions [21]. The production of a pseudoscalar Higgs via gluon 
gluon fusion is particularly simple because the leading order process is 2 ^ 1 and thus the kinematics 
is trivial. The case of scalar Higgs production has been studied in [21]; the structure of the result is the 
same for pseudoscalar case, the only difference being the values of the numerical coefficients computed 
in Tables 1 and 2. We write the rapidity distribution at the parton level as 



-{x,y;Tt,Tb) = 2^ -—{x,y,Tf,Tf>). 

fj'=t,b y 



dy 



(35) 
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rriA 


Cl{Tb,Tb) 


Cl{Tt,Tt) 


Cl{Tt,Tb) + Cl{Tb,Tt) 


100 


-1.9831 


2.8251 


-2.7106 


200 


-2.6344 


1.3629 


-3.8880 


300 


-3.0289 


0.35479 


-4.6358 


400 


-3.3123 


-0.37258 


-5.0118 


500 


-3.5334 


-0.46678 


-5.2424 


600 


-3.7147 


-0.55464 


-5.4629 


700 


-3.8681 


-0.63658 


-5.6627 


800 


-4.0014 


-0.71311 


-5.8438 


900 


-4.1190 


-0.78479 


-6.0089 


1000 


-4.2244 


-0.85211 


-6.1605 



Table 1: Results for the NLO coefficient ci, obtained by numerical evaluation of the integrals in Eq. (32). 



ruA 


C2iTb,Tb) 




C2{Tt,Tb) + C2{Tb,Tt) 


Cll{Tb,Tb) 


cii(rt,rt) 


CllijuTb) + Cu{Tb,Tt) 


100 


4.1527 


5.5035 


5.3683 


2.1777 


8.5243 


-0.9321 


200 


6.2393 


2.4237 


8.0566 


4.7678 


2.4575 


1.4378 


300 


7.7494 


1.5146 


10.1403 


6.7586 


0.78853 


4.0076 


400 


8.9476 


1.4378 


11.8815 


8.3770 


0.70446 


5.5309 


500 


9.9485 


1.4748 


13.1381 


9.7473 


0.47723 


6.1455 


600 


10.8118 


1.5279 


14.2203 


10.9403 


0.32497 


6.8722 


700 


11.5735 


1.5921 


15.1778 


11.9994 


0.22632 


7.6113 


800 


12.2566 


1.6637 


16.0402 


12.9539 


0.16702 


8.3396 


900 


12.8771 


1.7404 


16.8276 


13.8239 


0.13720 


9.0487 


1000 


13.4462 


1.8205 


17.5537 


14.6245 


0.12990 


9.7366 



Table 2: Results for the NNLO coefficients C2 and cn, obtained by numerical evaluation of the integrals 
in Eq. (33). 
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Figure 1: The NLO coefficient function in the quark-gluon channel Cqg in the case of just a top quark 
running in the loop, as a function of x. The green dashed curve is the effective theory result, while the 
sohd curves are full result [2] for nit = 172.0 GeV and niA = 100, 200, 300, 500 and 1000 GeV (from top 
to bottom) 
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The small-x limit of the different contributions are 



da, 



99 



Ay 



da, 



99 



dy 



da. 



dy 



{x,y;Tf,Tf 



TT 



6 ( y - ^Inxj + S (y + ^lnx 



Cl{Tf,Tf> 



+ 



^ [y- \ ^^^) + s (y + ^^^x 



In -C2{Tf,Tf>) 



+cii(r/,r// 



{x,y;Tf,Tf>) = 



1 + Sff 

2 



TT 



Cf 



6 ( y - ^ Inx j + 6 (y + ^lnx 



+ 



6 ( y ~ - Inxj + 6 (y + - \nx 



ci(t/,t/') 



1 C2(r/,T/0 



+cii(t) 



(36) 



4 Improvement of the NLO and NNLO cross-sections 

In the spirit of Ref. [13, 14] we construct an approximated cross-section by matching the coefficient 
functions computed in the effective theory cj^J'\x; oo) , which are vahd at large x, to the asymptotic 
behaviour at small x, computed keeping the mass of the fermions finite C^j ^Yneiiix 

{x;Tf). Wc concentrate 

first on the case of a top loop, which gives the dominant contribution at small tan p. We define 

n = l,2. (37) 

We start by considering the NLO case (n = 1), where the full result is already known. The spurious 
logarithmic growth of the effective theory at small x must be replaced by a constant behaviour. The 
matching point xq can be chosen at the intersection of the m/ ^ oo curve with the asymptotic constant. 
Subleading terms, which vanish at small-x, can be adjusted to ensure continuity. 

At the NNLO level (n = 2) the effective theory exhibits a cubic logarithmic behaviour, while the 
full theory has only a single logarithm. The approximated cross-section is constructed subtracting from 
the effective theory result all those terms which do not vanish at small x. However, at this order, the 
matching is more ambiguous because we only control the LLx behaviour, i.e. the coefficient of the 
single logarithm but we do not know the constant. Following [14] we choose to match the two curves at 
the point where the slope of the effective theory curve, on a logarithmic scale, approaches the correct 
asymptotic behaviour. The constant is then adjusted requiring continuity of the coefficient function. 



e(xo 



lim C^"'\x;(X)) 
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Figure 2: The NNLO coefficient function in the gluon-gluon channel Cgg in the case of just a top quark 
running in the loop, as a function of x. The green dashed curve is the effective theory result, while 
the solid curves are our matched results Eq. (37) with mt = 172.0 GeV and tua = 100, 200, 300 and 
500 GeV (from top to bottom). 



This method ensures that the approximated coefficient functions Eq. (37) and their first derivatives 
are continuous functions of x. As an example, in Fig. 2 we plot the NNLO coefficient function for 
the gluon-gluon subprocess, in the effective theory and in our approximation for different values of 
niA = 100, 200 and 500 GeV, in case of a top quark running in the loop. In case of small tan/3 we can 
neglect the contribution coming from 6-quark loops. In order to asses the numerical impact of finite 
top mass effect coming from the hard part of the partonic coefficient functions, we compute X-factors, 
defined as 

r^tf A a"'^^''^'"^{as]Xh,mA;mt) 

K {as;xh,mA;mt) = r , (38) 

a{as;xh,mA; oo) 

using MSTW08 NNLO parton distribution functions [22] for proton-proton collisions at y/s = 14 TeV. 
We obtain: 

= 1.000, mA = 100GeV, 
= 0.998, mA = 200GeV, 

= 0.985, mA = 500GeV, (39) 

The ET-factors clearly shows that at small tan /3, where the top quark contribution dominates, finite 
fermion mass effects are small in a wide range of the pseudoscalar mass, reaching a few percent only 
above the two-top threshold. 



12 



Our calculation of the LLx behaviour also includes the effect of a bottom quark running in the 
loop as well as the top-bottom interference. These contributions cannot be neglected at large tan f3. 
However, at first sight, it is not clear how to use this information in a NNLO matched calculation. The 
mass of the b quark is much smaller than niA and hence an effective theory where this contribution is 
integrated out is not justified. However, we notice that the large x behaviour of the coefficient function, 
being dictated by soft gluon effects, does not depend on the fermion mass. The effective theory result 
is exact in the x ^ 1 limit. This contribution dominates the gg A cross-section after convolution 
with the parton densities. It was indeed suggested that an effective theory where the bottom quark 
has been also integrated out could be, a posteriori, a reasonable approximation [23]. Therefore, we 
argue that a matching procedure like the one suggested for the top quark, should capture most of the 
contribution. However, we do expect the uncertainty to be larger as uncontrolled, mass-dependent 
terms at moderated x play a more relevant role. In order to asses the impact of finite fermion mass 
effect we compute K-factors, defined as the ratio of the matched inclusive cross-section to the effective 
theory result, for which top and bottom are both infinitely massive. For tan /3 = 10, we find: 



The K-factors show that, as expected, finite mass effects are more pronounced when we include the 
bottom in the loop. They are of the order of a few percent at low values of the pseudoscalar mass, 
reaching 5% for rriA = 500 GeV. However, we should stress that the matching procedure and the use 
of the effective theory at moderate x are less under control than in the small tan /3 case, where we only 
have to deal with the top quark. For this reason these iC-factor should be taken as an estimate. 

5 Conclusions and Outlook 

In this letter we have computed the leading small x singularities of the cross-section for pseudoscalar 
Higgs production in gluon gluon fusion to all orders in the strong coupling constant. We have also 
determined the small x behaviour of the Higgs rapidity distribution to the same accuracy. We have 
included the contributions of top and bottom quarks, together with their interference, showing that the 
cross-section only exhibits single high-energy logarithms, as opposed to the calculation in the effective 
theory where the quarks are integrated out. We have explicitly evaluated the coefficients of the leading 
logarithms up to NNLO. 

In the spirit of Ref. [13, 14], we have used our result to construct an approximation to the NNLO 
cross-section, by matching the effective theory result at large x to our asymptotic small x behaviour. In 
this way we can estimate the bulk of finite fermion mass corrections. At small tan /3, where the top quark 
contribution dominates, we have found that they are small, reaching 1-2% only for niA > 500 GeV. 
At large tan f3 the bottom quark contribution cannot be neglected. Including it in a similar way as we 
have done for the top, we have found an effect of a few percent. In order to refine this estimate, it 
would be important to match the small x behaviour to the asymptotic expansion in inverse powers of 
the fermion mass, of which the effective theory can be seen as the zeroth order term. 



^tb 



0.971, mA = 100GeV, 
0.964, mA = 200GeV, 
0.955, mA = 500GeV. 



(40) 
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